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Contribution to the Festschrift in honour of Professor Mario Novello. 

In this article we carefully distinguish the notion of bi-refringence (a polarization-dependent dou- 
bling in photon propagation speeds) from that of bi-metricity (where the two photon polarizations 
"see" two distinct metrics). We emphasise that these notions are logically distinct, though there 
are special symmetries in ordinary (3+l)-dimensional nonlinear electrodynamics which imply the 
stronger condition of bi-metricity. 

To illustrate this phenomenon we investigate a generalized version of (3+l)-dimensional nonlinear 
electrodynamics, which permits the inclusion of arbitrary inhomogeneities and background fields. 
[For example dielectrics (a la Gordon), conductors (a la Casimir), and gravitational fields (a la 
Landau-Lifshitz).] It is easy to demonstrate that the generalized theory is bi-refringent: In (3+1) 
dimensions the Fresnel equation, the relationship between frequency and wavenumber, is always 
quartic. It is somewhat harder to show that in some cases {e.g., ordinary nonlinear electrodynamics) 
the quartic factorizes into two quadratics thus providing a bi-metric theory. Sometimes the quartic 
is a perfect square, implying a single unique effective metric. We investigate the generality of this 
factorization process. 

PACS numbers: 



I. INTRODUCTION 



It is commonly assumed that electromagnetic phenom- 
ena (in sharp contrast with gravitational phenomena) can 
be described by a linear field theory even in strong fields. 
But there are many reasons to expect that this is not 
ultimately what happens. Born and Infeld 13] proposed 
that some specific non-linearities could appear at strong 
fields preventing the existence of arbitrary large values 
for the electric field surrounding a charged point parti- 
cle. For their part, Euler and Heisenberg [3] (see also 
Schwinger || ) argued that at high values of the electric 
field, the quantum creation and annihilation of electron- 
positron pairs would give rise to an effective non-linear 
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electrodynamic theory. This effective non-linear electro- 
dynamic theory can be further modified by externally 
driven alterations of the quantum vacuum owing to back- 
ground electromagnetic Q and gravitational fields [[|, 
or to geometrical boundary conditions (Casimir plates, 
etc.) H. For some additional important examples see 
the references in Q. Remarkably enough, in most cases 
these non-linear electrodynamic theories can be described 
by the use of appropriate (non-linear) media, with their 
associated refractive indices. 

Over the last few years Professor Mario Novello and 
co-authors have devoted considerable effort to investi- 
gating some of the consequences of having a non-linear 
electrodynamic theory. On the one hand, they have been 
interested in the possible role of non-linear electrodynam- 
ics in circumventing the singularity problem of general 
relativity. In particular, they have shown that in Euler- 
Heisenberg electrodynamics it is possible to have a cos- 
mology without an initial singularity Q. On the other 
hand, Professor Novello and co-workers have been inves- 
tigating what effects one should expect to see in the prop- 
agation of photons; now viewed as linear perturbations 
around a background electromagnetic configuration [p|. 
They were aware that in this case the causal propagation 
of photons is not controlled by the "physical" spacetime 
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metric (or others conformally related to it), but by an 
"effective" metric depending on the background electro- 
magnetic helds [[To). Professor Novello and his collab- 
orators realized that this opened up the possibility of 
building "geometric structures" in a manner analogous 
to, but different from, the realm of usual general rela- 
tivity. For arbitrary non-linear theories they have shown 
that black holes ]X1| , wormholes fll2| and even geometries 
with closed "timelike" curves |l3fl can be constructed. 
These "effective" geometries will only be felt by the pho- 
tons, while other matter fields will feel the usual gravita- 
tional spacetime metric. 



A very important point addressed by Professor Nov- 
ello and co-workers is the quite generic appearance of 
bi-refringence in non-linear electrodynamics 0. (Se e 
also the work of Plebanski fllO||, Dittrich and Gies p4[ , 
Schrodinger |l5| , and Boillat [p6|.) The two polarization 
states of the photon propagate differently. The present 
authors have similarly been confronted with the question 
of bi-refringence (or more generally multi-refringence) in 
general systems of second order partial differential equa- 
tions (PDEs); this investigation bein g m otivated by an 
abstract approach to analog relativity [17) . (The last few 
years have seen a proliferation of analog models of/for 
general relativity in the literature. See |l8[ for an ex- 
tensive reference list.) In a series of papers |l7], |l9| we 
have shown that the crucial issue in building an analog 
model of general relativity is the linearization of non- 
linear field theories around some background solution. 
In the present article we will apply the general analy- 
sis and language of |l7| to a generalization of non-linear 
electrodynamics. We will show that the existence of bi- 
refringence is quite easily established, but that the step 
to bi-metricity (the existence of two different effective 
metrics controlling the propagation of each photon po- 
larization) requires special conditions that are satisfied 
by electrodynamics in (3+1) dimensions. 



We would also like to point out that nonlinear electro- 
dynamics (in particular Born-Infeld theory) has in recent 
years seen a marked resurgence of interest with the ad- 
vent of the notion of D-brane (see for example Polchin- 
ski I^TJ). Many physicists feel that D-branes will be a 
crucial ingredient in any final formulation of M/String 
theory. It happens that the motion of a D-brane in the 
bulk spacetime is controlled by a Born-Infeld type ac- 
tion |2l[] . This implies that while closed strings propagate 
following the bulk spacetime metric, open strings (whose 
end points are attached to D-branes) follow an effective 
metric derived from the Born-Infeld Lagrangian |25| . 



Finally, we point out that much of the formalism de- 
veloped below owes a great debt to related work (by SL, 
Sebastiano Sonego, and MV) on photon propagation at 
oblique angles in the Casimir vacuum |23j. 



II. GENERALIZED NONLINEAR 
ELECTRODYNAMICS 



Consider a general class of Lagrangians of the form 

^effective = £ (F^(x), B(x)) . (1) 

Here denotes the electromagnetic field strength; and 
it is assumed that derivatives of this field strength do not 
occur in the Lagrangian. In terms of the vector potential 



(2) 



In addition B{x) denotes a generic class of external non- 
dynamical background fields. These could represent, for 
example, a refractive index, the 4- velocity of a dielectric, 
the location and 4-velocity of Casimir plates or other 
conductors, assorted inhomogeneities and/or boundary 
conditions, an external gravitational field, etc. If these 
background fields are all set to their trivial position- 
independent values then the system reduces to ordinary 
nonlinear electrodynamics in which the Lagrangian de- 
pends only on the two independent Lorentz invariants 
that can be constructed from the field strength tensor 
(for example, Born-Infeld or Euler-Heisenberg electro- 
dynamics). 

The complete equations of motion for nonlinear elec- 
trodynamics consist of the Bianchi identity, 



plus the dynamical equation 



0, 



(3) 



(4) 



We now adopt a linearization procedure: Split the elec- 
tromagnetic field into an internal (possibly dynamical) 
background field plus a propagating photon 



P ^background _r_ ^photon 



(5) 



Then, assuming the background satisfies the equations of 
motion and retaining only linear terms in the propagating 
photon, we have 



(/ photon )[^a] = 0, 



(6) 



and 



d 2 C 



dF^ dF aP 



photon j = Q _ (7) 

background / 



On defining 



d 2 C 



dF^ dF a p 



(8) 



background 



equation (^) can be rewritten in the somewhat more com- 
pact form 



0, 



(9) 
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Note that the tensor Vt^"" is symmetric with respect to 
exchange of the pairs of indices \xv and a/3, and antisym- 
metric with respect to exchange of indices within each 
pair. That is: Sl^ va P has most of the key symmetries of 
the Riemann tensor. If one wishes to work directly at 
the level of the linearized Lagrangian one has: 

f jJ_Q/ii/Q/3 jphoton /•photon 

^-linearized — g Jfif ' af} ' \ lyJ ) 

The key observation is that this linearized Lagrangian 
generically leads to birefringence. 

We now apply a restricted form of the eikonal approx- 
imation by introducing a slowly varying amplitude f^ v 
and a rapidly varying phase <j>: 

f$° ton = e 1 *. (11) 

The 4-wavevector is then defined as k^ — d^cf). This 
approximation is similar to, but not quite identical with, 
the usual eikonal approximation. This is because one 
assumes that 4> varies on scales much smaller than those 
of the background, while, on the other hand, use of the 
Lagrangian (Q) also implies that the components of k are 
much smaller than the values fixed by the electron mass. 
(This is the so-called soft-photon regime). Under these 
hypotheses, 

SV tavf> k a f vp = 0. (12) 

But in general the internal dynamical background field 
is itself subject to quantum fluctuations, and to take this 
into account the coefficients of this equation are to be 
identified with the expectation value of the corresponding 
quantum operators in the background state \tp): 

(il>\W a »^)k a f vP = Q. (13) 

In taking this expectation value we are using the fact that 
the fluctuations in the internal dynamical background 
fields are influenced by the external non-dynamical back- 
ground fields B(x). (For example, in the case of the 
Casimir geometry, by the distance between the plates.) 
In the spirit of the restricted eikonal approximation there 
is a separation of scales between the internal dynamical 
background fluctuations and the propagating photon. 

The Bianchi identity (@) constrains to be of the 
form 

f^u — kfi a v k v dp, (14) 

where we have introduced the linearized gauge potential 
a for the propagating field. Inserting ( p^ ) into ( |l3| ) we 
find 

(fl>\Sl" m, ^\i>)k a kpa t/ = 0. (15) 

Note that this last equation implies that any completely 
antisymmetric part of (^ val3 ) can be discarded without 
affecting the equations of motion. 



Of course, this entire discussion could alternatively be 
rephrased in terms of Hadamard's theory of the propaga- 
tion of weak discontinuities pj| , the formalism preferred 
by Professor Novello |J . An identical equation (relating 
the polarization and the wavevector) is encountered. 

We emphasise that the discussion in this section has 
(so far) been completely independent of the dimension- 
ality of spacetime. If we now ask how many independent 
components arise in 

{^\n^ a0 \ip) = (^(MM) ip) (16) 

we find [in (d+1) dimensions] 

ri(rf+l)(ri 2 + rf + 2) 
8 

In particular in (3+1) dimensions this quantity has 21 
independent components. One of these components can 
be taken to be the coefficient of the Levi-Civita tensor, 
and so does not affect the equations of motion ( [TBI) . An- 
other component can be interpreted as the overall scale 
of f2, which again does not affect the equations of motion. 
So the number of useful independent components in Q is 
19. In contrast, two light cones only specify 2 x 9 = 18 
components. (Two metrics would specify 2 x 10 = 20 
components, but in this paper we are only looking at the 
null cones.) It is ultimately this close relationship which 
makes (3+1) dimensions so special (and tricky). 

III. FRESNEL EQUATION 

Equation ( |l5|) represents a condition for a as a function 
of k — it constrains a to be an eigenvector, with zero 
eigenvalue, of the fc-dependent matrix 

A^(k) = k a k p . (18) 

Any non-zero solution corresponds to a physically pos- 
sible field polarization, that can be identified by a unit 
polarization vector e (provided a is not a null vector — 
a possibility that can always be avoided by a suitable 
gauge choice). 

A necessary and sufficient condition for the eigen- 
value problem A^ u a u — to have non-zero solutions is 
det(A tl ' / ) — 0; however, this gives us no information at 
all. Indeed, any a parallel to k is always a non-zero solu- 
tion, so the condition det (A^) = is actually an iden- 
tity. On the other hand, a \\ k is merely an unphysical 
gauge mode that corresponds to f^ v = by (|l4j), so we 
need to find other, physically meaningful, solutions of the 
eigenvalue problem. 

To this end, we exploit gauge invariance under a — > 
a+X k and fix a gauge, thus removing the spurious modes. 
For the current subsection, it is particularly convenient 
to adopt the temporal gauge eto = 0. Then we can de- 
fine a polarization vector e M = a^j (a^a^) 1 ^ 2 , and the 
eigenvalue problem A^ v e„ = splits into the equation 

A 0i e t = 0, (19) 
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plus the reduced eigenvalue problem 



A^e j = 0. 



(20) 



The latter admits a nontrivial solution only if 

det (A ij ) = 0. (21) 

The condition ( ^l|) plays the same role as the Fresnel 
equation in crystal optics [ p4| — it is a scalar equation 
for k and thus gives the dispersion relation for light prop- 
agating in our "medium" . 

The spatial components of the matrix are 

A 2j = lo 2 (4>\n t0l0 \i/j) + ojk m (rjj\n i03m + n imj '°|i/'> 
+k m k n (^\n im ^\i>), (22) 

where lu — fco- Let us define the unit vector k = k/\k\. 
Then the components of A %3 in a basis with one axis 
directed along k are 

A 13 kj = lu 2 (^|rr 0j0 |V>) kj+uk m (^|O <mi0 |^> k 

= lu (u (^p\n i0: > \^) kj + k m (ip\n im: > \^) kj 

= luV\ (23) 
In particular 

A ij k t k 3 = lu 2 (ip\Cl iOj0 \ip) kk = lo 2 S. (24) 

If we now specialize for dcfinitcncss to (3+1) dimensions, 
then the matrix A 13 has the following structure 



lu 2 S luV j 
lu V 1 T IJ 



(25) 



where / and J label the two directions orthogonal to k 
in the sense of vector-space duality. (The V 1 are linear 
in the 4-wavenumber k, while the T IJ are quadratic in 
the 4-wavenumber k.) Evaluating the determinant by 
expanding in the first row or column, it is easy to see 
that every term will contain at least two factors of lo, 
which establishes 



det (A ij ) = lo 2 P 4 (k), 



(26) 



where V±(k) is a homogeneous fourth-order polynomial 
in the 4-wavenumber k^ = (u>,\k\ fcj). While the deter- 
minant is itself a sextic, the physically interesting part is 
given by the quartic V±. In fact, by the rules for parti- 
tioning determinants 



det (A ij ) =lo 2 S det 



rpIJ _ 



V 1 V J ~\ 



so that V± is effectively a 2 x 2 determinant 

V 1 V J ~ 



Vi{k) = S det 



T 



u 



S 



(27) 



(28) 



where the elements of the 2x2 matrix are themselves 
quadratic in the 4-wavenumber. This allows us to re- 
phrase the current analysis in the language of our more 
general "normal modes" analysis jtTj by introducing a 
matrix f^' IJ (extremely similar to but not identical with 
the related quantity introduced in that article) and writ- 



ing 



Vi(k) = det [/""'" fc M K] 



(29) 



Here the two polarization states take on the role of (dual) 
"field indices" as discussed in reference jD]]. The deter- 
minant is to be taken on the I J indices. 

The upshot of this analysis is that in the most general 
case there appear to be four dispersion relations, corre- 
sponding the four roots of the quartic equation 



Vi{k) = 0. 



(30) 



If we write k = (lu, k) = (lo, \k\ k) then two of these roots 
correspond to propagation in the +k direction, while the 
other two correspond to propagation in the —k direction. 

Different polarization states are represented by linearly 
independent solutions of the eigenvalue problem (20), un- 



der the condition (21 ). Thus, the space of polarizations is 
exactly two-dimensional. Since equation (|2l]) gives rise to 
two dispersion relations, the polarization states actually 
satisfy two (in general, different) eigenvalue equations, 



e W=0, 



(31) 



where r — 1,2 labels the dispersion relations and A^ v is 
obtained from A^ v by imposing the corresponding con- 
dition on k as derived from equation (|30|). Indeed, sup- 
pose you pick a specific 3-direction k and have by some 
means determined two independent polarization states 

(r) 

e v , which are implicitly functions of k and the corre- 
sponding solutions k of Vi(k) = 0, then one can construct 
a pair of two-index matrices 



3& = {^W e£ ^ 



Ar) >) 



(32) 



Although these quantities are matrices with the correct 
index structure to be interpreted as "effective metrics" 

they are in the general case implicitly functions of lu, \k\, 
and the direction k, and so these quantities cannot be 
viewed as spacetime metrics. 

It is only in some special cases (e.g., ordinary nonlin- 
ear electrodynamics) that the polynomial V&(ki) factor- 
izes into two quadratic forms, 

V A (k) = (gft k^ k v ) (Sg k a kp) , (33) 

in which case we obtain two second-order dispersion re- 
lations: 



and 



^(2) k^ky — 0, 



(34) 
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with momentum-independent matrices Qua. We can now 
interpret these two matrices Q as the two (inverse) ef- 
fective metrics of a bi-metric theory. (More precisely 
they are representative elements of two conformal classes 
of inverse metrics, since multiplication by an arbitrary 
position-dependent scalar will not modify the dispersion 
relations.) 

In very special cases not only does (E53l) hold, but also 



f(i) 



G(2)- (That is, the fourth-order polynomial Vi(k) 



is a perfect square.) In this case one ends up with a single 
quadratic dispersion relation of the familiar form 



k v 



0. 



(35) 



where is some symmetric tensor, which we shall call 
the (inverse) effective metric (again defined only up to 
an arbitrary conformal factor). 

We now wish to investigate the conditions under which 
these factorization (bi-metric) and uniqueness properties 
hold. Since the entire formalism ultimately derives from 
the tensor r^ aiy/3 , we shall look for suitable algebraic con- 
straints on this tensor. 



IV. SINGLE EFFECTIVE METRIC 



In the case of a single effective metric we have 



g" v k K = o. 



(36) 



It should be clear from our previous discussion that a nec- 
essary condition for this to happen is the absence of bire- 
fringence. We see that the wave vector is now null with 
respect to this (unique) "effective metric" , which 
therefore defines an effective geometry for the propaga- 
tion of light. 

We warn the reader that even when a unique (inverse) 
effective metric Q^ v is defined, we always raise and lower 
indices using the flat Minkowski metric tj^u, or in the 
presence of a gravitational field the physical spacetime 
metric g^. The effective metric itself, denoted g^, is 
the matrix inverse of Q^ v . Because of the way indices are 
raised and lowered using the physical metric you cannot 
use index placement to distinguish Q = g -1 from g. 

This single effective metric situation implies that (up 
to possibly a piece proportional to the Levi Civita ten- 
sor) the tensor (■0|f2 A ' O!1 ^ 3 |?/>) must be algebraically con- 
structive solely in terms of G^ v . In view of the symme- 
tries of fl tJ - c " , f 3 we know, without need for detailed calcu- 
lation, that it must be of the form 



{ij}\Qf iavP \4i) = ^(g^ v g aP - g va ) + $ 



(37) 



for some quantities W and By appealing to the confor- 
mal invariance of the null cones we can always absorb a 
factor of 2^/|^| into the inverse metric g and so rewrite 
this as 



.fcg/iautf ( 38 ) 



Conversely, if f^ Qly/3 is of the form then the matrix 
©is 



{g^{g af) k a k p ) -(g" a k a ) (<r%)}, 



(39) 



and the photon propagation equation, A^ v a v = 0, be- 
comes 

(g aP k a k p )g^ a v - (g al3 a a k p ) Q» v k v = 0. (40) 

This equation is obviously satisfied by the uninteresting 
gauge modes a || k, with no constraints on k. Solutions 
corresponding to a non-vanishing f^" exist only if the 
coefficient of g^ v a v is zero, i.e., if ( |35| ) holds. Thus, the 
two polarization states propagate with the same disper- 
sion relation (p5), and there is no birefringence. 

Substituting ( 35| ) back into the propagation equation 
( [f"o| ) we find another relationship typical of this case, 



g^ K a. 



0. 



(41) 



Formally, the above equation looks like a gauge condi- 
tion. This might seem puzzling, because nowhere in the 
present subsection have we fixed a gauge. In fact, ( |4l] ) 
is a consequence of the dynamical equation A* 1 " a v = 0, 
when the "on-shell" condition ( p5| ) is satisfied, and it does 
not imply any gauge fixing. 

It is also interesting to notice that there is now a self- 
consistency or "bootstrap" condition, 



(42) 



We stress that these relations depend only on the as- 
sumed existence of a single unique effective metric g^^ — 
they do not make any reference to other specifics of the 
external background fields B(x) or the quantum state. 

Finally we point out that in this mono-refringent case 
the linearized Lagrangian reduces to 



linearized 



qPl 1 y phot on qvol yphoton 



~* ^ f$ oton C oton - ( 43 ) 



With hindsight, this is exactly what we should have ex- 
pected. If we now use this Lagrangian formulation to 
demand positivity of energy [a feature missing from the 
purely kinematical analysis based on equation ([Hf)] then 
we should set ± — > +1. Note that we have also used the 
conformal invariance of the null cones to normalize g in 
the conventional manner. Finally the $ term is simply 
the well-known Pontryagin index. 



V. ORDINARY NONLINEAR 
ELECTRODYNAMICS 

In order to see how to develop a general ansatz that 
leads to bi-metricity (perhaps not the most general 
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ansatz) it is useful to consider the explicit form of the 
tensor f2 M „ Qj a for ordinary nonlinear electrodynamics: 

£ NLE = £(T,G). (44) 

Here we have adopted the now common variables ]l4| 

1 „ „„„ 1 



T = j-Fjui/ F**" = -[B — E ) , (45) 



Q = ~F IJV *F' tv = -E-8. 



(46) 



For such Lagrangians 



1 



+ - {d g C) e^ af3 



+F> iV F al3 (d%£) + *F" V *F af} (d 2 g C) 

+ (pa* *p a p + *p»» pa ^ d:Fg£ _ ^ 

As soon as one inserts this tensor into the photon equa- 
tion of motion (|lj), the completely antisymmetric part 
proportional to the Levi-Civita tensor drops out, because 
of the Bianchi identity (^). The remaining pieces repro- 
duce the photon equation of motion in the perhaps more 
usual form considered by Dittrich and Gies jl4|, or Nov- 
ello and co-workers 

Unless one has a specific need to perform calculations 
to orders higher than 0(a 2 ), it is often sufficient to con- 
sider the Euler-Heisenberg Lagrangian which, in the 
T-Q formalism adopted above, takes the form 



£eh = ~ F + Cl F 2 + c 2 g 2 , 
Air 



with 



Cl 



or 



la 1 



90n 2 mi 



(-2 



360tt 2 to4 • 



(48) 



(49) 



The terms proportional to !F 2 and g 2 of this Lagrangian 
are quartic in the field, and describe the low-energy limit 
of the box diagram in QED, when four photons couple 
to a single virtual electron loop. Thus, the Lagrangian 
( [48| ) is only accurate to order a 2 , and it is meaningless 
to retain higher order terms within this model. For the 
Euler-Heisenberg Lagrangian, the tensor Q^ va P is 

1D7T 2 1 ' 



^2 g r fiua/3 



1 pliv pa(3 , £2 * pfif * pap 

2 2 



(50) 



Again, when one inserts this tensor into the photon equa- 
tion of motion Jl^), the completely antisymmetric part 
proportional to the Levi-Civita tensor drops out. 
Suppose we now adopt the "rotated" quantity 



K = cos 6 F + sin 6 *F, 



(51) 



so that 



sin 6 F + cos 6 *F. 



(52) 



Then for both (ordinary) NLE and Euler-Heisenberg 
electrodynamics we can put fi into the form 



+<j> 

+a K" v K a/3 



p/3 ^VOL 



b*K" u *K al3 . (53) 

The matrix is then [suppressing explicit indices] 

A = ^{rj{k,k) 7] - (rjk) ® (nk)} 

+a{Kk)®{Kk) + b{ i 'Kk)®( k Kk). (54) 

Note that because both K and *K are antisymmetric 
K(k, k) and *K(k, k) = 0; therefore A k = as ex- 
pected. [We have defined r](k, k) = rj^ v k^ k v , K(k, k) = 
K^ v k v , etc] Now consider det'(A), the "reduced" 
determinant in the three directions orthogonal to k. We 
adopt this particular reduced determinant as a techni- 
cally convenient alternative to using the temporal gauge. 
To be precise we define 



det'(A) = 



d[det(A + e/)] 



de 



(55) 



c=0 



which implies 



1 



det'A = ^{[Tr(A)] 3 + 2[Tr(yL 3 )] - 3[Tr(A)][Tr(A 2 )]}. 



(56) 



Now use this formula, or the fact that in 3 dimensions 



det 3 (AJ 



u®u- 

i2 



v <g> v) 



(57) 



= A {A 2 + \(u 2 +v 2 ) + [u 2 v 2 ~{u-v) 2 ]} 
(useful when k is timelike with respect to 77), to deduce 

det'(A) = -*I(Jfc,fc) |* 2 r/(fc,/c) 2 

n(k, k) [a j(Kk, Kk) + b j(*Kk, *Kk)} 
+ab[-f{Kk, Kk) j(*Kk, *Kk) 

-j(Kk,*Kk) 2 ]y (58) 

[The presence of the — sign arises from the indefinite na- 
ture of the metric 77. We have also verified the above 
formulae via explicit evaluation of the determinants us- 
ing Maple. Note that we now distinguish between the 
(inverse) contravariant Minkowski metric r\ and the co- 
variant Minkowski metric 7 = vj~ . The wavevector 
k is always taken to be covariant while both K and 
*K are assumed doubly contravariant. We have defined 
j(Kk,Kk) = 7 M „ {KkY {Kkf, etc.] 

Discarding the uninteresting factor of I(k, k) = k T k 
(it corresponds to the factor lo 2 encountered when we 
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worked in temporal gauge) we identify 

V 4 (k) = * 2 7](fc,fc) 2 

r)(k, k) [a j(Kk, Kk) + b 1 {*Kk, *Kk)} 
+ab[j(Kk, Kk) j(*Kk, *Kk) 

-^(Kk^Kk) 2 }. (59) 

This is clearly a quartic in fc, and the "miracle" of (ordi- 
nary) nonlinear electrodynamics is that it factorizes into 
two quadratics. To establish this factorization we use 
some very special properties of (3+1) dimensions: 

Tm^wcII = F 7 F - - tr (F7F7) n (60) 

= * F 7 *p _ tr (*F 7 *F<y) 77 (61) 

= *F 7 *F + - tr (F7F7) 77 (62) 

and 

F^*F = Qr/ = *F^F. (63) 

In terms of the "rotated variables" the Maxwell stress 
energy tensor is given by 

TMaxwoii = K 7 K - - tr (R^Kf) n (64) 

= * K 7 *K - - tr (*K 7*Jf7) ry (65) 

= 7 *if + - tr (KjKj) r). (66) 

Additionally 

K 1 *K = g e ri = *K 7 if, (67) 

where 

J" e = cos(26») + sin(20) £ (68) 

Q B = -sin(20) .F + cos(26>) £. (69) 

From these relations we can deduce 

j(Kk,Kk) = -T Maxwea (k, k) - T e r)(k,k), (70) 
7 (*Kk,*Kk) = -T Maxwe n(fc, fe) + ^ e tj(fc, fc), (71) 
j{Kk,*Kk) = -g e v(k,k). (72) 

When substituted into V 4 this implies (dropping the ex- 
plicit "Maxwell" subscript) 

V 4 = a r](k, k) 2 + ai T)(k, k) T{k, k) + a 2 T{k, k) 2 , (73) 

for suitable ao, 01, a2- Indeed 

a = * 2 + *(a-6)J' e -a&(V + S<> 2 ) , (74) 
01 = #(a + 6), (75) 
a2 = a&. (76) 



The quartic will factorize provided 

ao + aix + CL2X 2 — (77) 

has a solution in the real numbers. Fortunately the dis- 
criminant is easily seen to be a sum of squares and so is 
always positive. 

a\-Aa a 2 = [*(o - 6) + 2 afejg] 2 + [2 ab Q d f > 0. (78) 

This is now enough to guarantee that V 4 factorizes, in- 
deed 

V 4 = a [ v (k, k) + h T(k, k)} [ V {k,k) + b 2 T{k,k)} (79) 

for suitable ao, 61, b 2 . These coefficients will be functions 
of ^ , a, b, Tq and Qg whose precise form is not needed 
for the point we are currently making: As long as the 
Lagrangian is only a function of the two invariants T 
and Q, then the theory is not just birefringent, it is truly 
bi-metric. 



VI. A GENERAL BI-METRIC ANSATZ 

Based on the above we can now guess a general ansatz 
for fi^"^ that always leads to bi-metric propagation (we 
do not guarantee that this is the most general ansatz). 
Let 

= * {f° g vp - g^ g va ) 

+a K af3 +b*K» v *K af3 , (80) 

where g^ v is now any symmetric matrix of Lorentzian sig- 
nature and K^ v is any anti-symmetric matrix — in par- 
ticular K does not necessarily have anything to do with 
the electromagnetic field. All the analysis in the previous 
section can now be converted into purely algebraic state- 
ments about g and K . For example T no longer has the 
interpretation of being the Maxwell stress-energy tensor, 
it is simply an algebraic matrix defined by 

T = Kg- X K - -tx{Kg- l Kg~ l )g. (81) 

Reinterpreting everything in this purely algebraic man- 
ner, and repeating the analysis of the previous section, 
we see that ( p0| ) leads to bi-metric propagation with the 
two light cones being given by linear combinations of g 
and T. 

Though rather general, this is not likely to be the most 
general bi-metric ansatz. To see this note that (|^) above 
appears to contain 10(g) + 6(K) + 4(\E r , $, a, b) = 20 free 
parameters. But 4" can be absorbed by redefining g, 
while $ does not affect the equations of motion, a can 
be absorbed by redefining K, and an overall scale factor- 
does not affect the equations of motion. This leaves 16 
physically interesting free parameters in ( |80| ) to be com- 
pared with 2x9 = 18 free parameters encoded in a pair 
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of light cones. We have not as yet been able to deduce the 
most general form of VL^ va @ compatible with bi-metricity. 

In terms of the linearized Lagrangian our bi-metric 
ansatz ( |so|) corresponds to 

^linearized ~Z ^ (/photon 9 /photon 9 ) 

+ \® tT (/photon g' 1 * /photon 
+ \ atr ( K 9~* /photon 

+ 1 b trCKg- 1 /photon^ 1 ) 2 - (82) 

VII. BIREFRINGENCE WITHOUT 
BI-METRICITY 

To wrap up, can we now give a simple explicit example 
of a model that is birefringent without being bi-metric? 
Suppose we have a pair of two-forms Ji and J 2 with 
Ji gt*J 2 - Then consider (as a particular example) 

+a Jf Jf^ + b J£ v Jf. (83) 

In this case, because you are now not satisfying the spe- 
cial algebraic constraints of the previous section, there is 
no reason for the determinant to factorize. Indeed 

V 4 (k) = ^ 2 n(k,kf 

n(k, k) [a 7(Jifc, Jik) + b 7 ( J 2 k, J 2 k)} 
+ab[j(Jik, Jik) 7(J 2 fc, J 2 k) 

-~f(Jik, J 2 k) 2 }. (84) 

In terms of the linearized Lagrangian (suppressing the 
inverse Minkowski metric 7 = as being understood) 

^linearized ~^ ^ ^ (/photon) ~h ^ (/photon /photon) 

+ \ {atr(J X /photon) 2 +&tr(J 2 /photon) 2 } • 

(85) 

The key point here is that the linearized Lagrangian is not 
a simply a function of the invariants T and Q. In going 
to generalized nonlinear electrodynamics we have permit- 
ted additional structure in the form of external fields and 
boundary conditions. Unless these external fields satisfy 
rather particular conditions (e.g., in the present example, 



Ji oc * J 2 ) there is no reason to believe the Fresnel de- 
terminant factorizes, and no reason to expect a bi-metric 
theory. 

Indeed, if one picks a tensor f2 "at random" and ex- 
plicitly evaluates Pi(k) (using a symbolic program such 
as Maple) one rapidly concludes that arranging factor- 
ization (and so bi-mctricity) is not an easy task. Bi- 
metricity is not generic in the set of all birefringent the- 
ories. 

VIII. CONCLUSIONS 

In this Festschrift article we have discussed, on 
quite general grounds, the phenomenon of bi-refringence 
and bi-metricity in [generalized] (3+l)-dimensional non- 
linear electrodynamics. Our treatment encompass any 
non-linear electrodynamic theory in interaction with an 
arbitrary number of external non-dynamical fields char- 
acterizing a general medium (for example, a flowing di- 
electric, external gravitational fields, moving Casimir 
plates, etc.). We have seen that the phenomenon of 
bi-refringence is both generic and easily established in 
these theories. In (3+1) dimensions, the Fresnel equa- 
tion is quartic and in special cases (e.g., ordinary nonlin- 
ear electrodynamics) factorizes into two quadratics (i.e., 
two metrics), typically different from each other. In more 
specialized situations these two metrics can be identical, 
leaving no opportunity for bi-refringence. 

If we consider generalized nonlinear electrodynam- 
ics, then because of the presence of additional back- 
ground fields, the close link between bi-refringence and 
bi-metricity can be broken — in such situations the Fres- 
nel equation is intrinsically quartic and to naturally de- 
scribe the geometry one would need to go beyond the 
notion of Lorentzian geometry and instead introduce the 
notion of a pseudo-Finsler geometry as described in [jl7| . 

In closing we emphasise that the use of nonlinear ex- 
tensions to electrodynamics is currently becoming ubiq- 
uitous. The implied notions of birefringence, bi-metricity, 
and pseudo-Finsler geometries will doubtless continue to 
attract considerable attention. Professor Novello's work 
on effective geometries and birefringence will continue to 
have important repercussions down the road. 
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